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Abstract. The speciality index, which has been mainly used in Numerical 
Relativity for studying gravitational waves phenomena as an indicator of the 
special or non-special Petrov type character of a spacetime, is applied here in the 
context of Mixmaster cosmology, using the Belinski-Khalatnikov-Lifshitz map. 
Possible applications for the associated chaotic dynamics are discussed. 



PACS numbers: 04.20Cv 

1. Introduction 

The speciality index (SI) ^ is an invariant and dimensionless indicator of tlie special 
or non-special Petrov algebraic character of a given spacetime. This quantity has 
been mostly used in Numerical Relativity, for instance to study a black hole which 
is radiating gravitational waves, or black holes merging E]. In this Letter we 
use the SI in the context of Cosmology, to study the Belinski-Khalatnikov-Lifshitz 
(BKL) map O^IIS], defining in this way an invariant indicator of the Petrov type 
"fluctuations" close to the singularity of the Mixmaster solution. Possible applications 
in the context of the associated chaotic dynamics are discussed. 

2. Petrov classification: speciality index. 

The algebraic properties of curvature are a very useful tool to obtain powerful insights 
into the character of a given spacetime metric. In particular those of the Weyl tensor, 
which is the trace free part of the Riemann tensor and in vacuum coincides with 
it, play a central role in Einstein's General Relativity theory. The Penrose-Debever 
equation l[uCpjqr[slt]^'^l'^ — states the existence of four distinct null eigenvectors for 
the most general spacetime: these are known as ^'principal null directions" (PND) [S] . 
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If all the PND result distinct one has the algebraically general case (Type I). When 
some of them coincide, this gives rise to the algebraically special case summarized as 
follows: Type II (one pair of PND coincides), Type D (two pairs of PND coincide). 
Type III (three PND coincide). Type N(all four PND coincide) and Type O (no PND, 
because of conformal flatness). 

Defining the complex tensor Cabcd = Cabcd — i*Cabcd, one can introduce the two 
complex curvature invariants 

J ^abcd J f^cd /^mnab fO ^\ 

-< — 'TT^'^abcd'^ J 'J — oo/l ^".bcd'^ mn'^ ■ 

These can be used to define the speciality index jT] 

S = 27J^/I^ (2.2) 

which marks, in an invariant way, the transition from certain algebraically special 
solutions (5 = 1) and the general Petrov type / {S ^ 1) ;:6^. We point out that for 
some spacetimes this quantity might be not well defined ,2. because of the possible 
vanishing of / and/or J, although for the vacuum Kasner spacetime 0, as shown in 
the following, this is not the case. This solution of Einstein equations is given by 

ds^ = d<2 - t^P'dx^ - t^P^dy^ - t^P'dz^ , (2.3) 

where 

Pi+P2+P3=pi+Pl+Pl^l (2.4) 

and the indices can take values in the interval [— ^ , 1] only. The Kasner metric admits 
two special subcases when two of the pi indices are equal: it then follows from H2.4|l 
that either pi = P2 = 0, = 1 (and permutations) and the spacetime is flat in 
this case, or pi ~ —1/3, p2 ~ Ps = 2/3 (and permutations) and one has the Kasner 
locally rotational symmetric type D solution, with a spindle-like singularity [2||Sj. For 
other choices of the parameter the Kasner spacetime is of Petrov type I. A simple 
calculation, using the constraints listed above, shows that the SI for the vacuum 
Kasner spacetime is 

27 27 

S = -~^PiP2P3 = -jpI{1-P3) ■ (2.5) 

By direct inspection this quantity ranges from the value 5=1 (the type D case) to 
S = Q (the flat case) and is well defined with continuity for any value of the parameter 
P3- 

3. Applications to BKL dynamics 

The class of algebraically general Kasner spacetimes contains only algebraically special 
subcases of either type D and type O, and its SI has the simple time independent 
form (|2.5|) . This allows us to use the latter to describe the Mixmaster dynamics 
as approximated by BKL Kasner epochs [31 0] |51 [7]. By using the standard 
parametrization: 

U M + 1 uiu +1) 

Pi = J- — T' P^ = ~^ TT' P3 = — T' we[l,+oo) (3.1) 

satisfying the ordering 

-^<Pi<0, 0<P2<| ^<P3<1, (3.2) 
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with the sequence of Kasner epochs given by the rule (Gauss map): 

Un+l =Un-l if 2 < U„ < OO , 

Un+l = ^—r if 1 < < 2 , 

Un - 1 

we obtain from (|2.5() an "n" -dependent speciality index labelled by each epoch 
27 uliun + 1)2 



(3.3) 
(3.4) 

4 (ul + un + ir- ^^-^^ 

Although irrational initial values ui only should be considered j7|, in numerical 
simulations this requirement results clearly an abstraction, because truncated rational 
numbers only can be handled. Using as an example the sequence given by Berger |1U| 
which starts "close" to the flat spacetime configuration with ui — 7.2328, we get the 
speciality index evolution shown in Figure JQ). 




Figure 1. Speciality Index for the vacuum BKL map with ui = 7.2328 in 
terms of epochs n. Interpolation for visual clarification only. The flattening 
on the right is due to the map sensibility. 

We point out that as soon as the system gets close to the type D configuration, 
it moves rapidly towards to the type O region, and then gently evolves until it gets 
close to the type D case again, and so on. The flattening of the trend in Figure 
is due to the well known sensibility of the Gauss map which can generate very long 
periods of monotonic behavior until oscillations start again. 

Pictorially, using the representation of the Mixmaster as the motion on a 
contracting triangular potential well in the time direction moving away from the initial 
singularity, the motion close to the type D case corresponds to an almost perpendicular 
bounce on the middle of the side of the triangle, with the incoming free-motion phase 
before the bounce corresponding to the flat space Kasner indices and the outgoing 
free-motion phase after the bounce corresponding to the Kasner indices. 

The bounce is equivalent to a transition from the one set of Kasner indices to 
the other, for the asymptotic behavior away from the straight wall but still far from 
its time spent at the opposite end in the "channel" corner of the potential where the 
exact Taub solution originates and then returns, and space curvature effects remain 
important since the system point is always close to the potential walls of the channel. 
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We point out that the BKL parametrization has cancelled any information 
concerning the specific direction in which the motion is happening, leaving in our 
case an invariant dynamics in an "abstract Petrov space" . 

4. Conclusions 

We have introduced the speciality index in Mixmaster BKL dynamics, in analogy with 
its use in the numerical treatment of gravitational wave sources. Because of its gauge 
invariant nature, time independence and adimcnsionality, the Kasner SI and its derived 
BKL version can be used in the sophisticated numerical simulations of Mixmaster, to 
get useful invariant information concerning chaos ^2 ^] ■ In particular it would be 
useful to study the probability associated on the various regions of the segment [0, 1] 
in which the SI ranges during its evolution. As pointed out in a recent review on the 
subject J3] (see also references therein), "a remaining open question is how closely 
an actual Mixmaster evolution is approximated by a single BKL sequence." A direct 
comparison of the BKL SI (which approximates the Mixmaster dynamics) with the 
corresponding exact Bianchi IX one will be the appropriate way for approaching the 
problem numerically, to give an answer to this question. 
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